Break-down of the density-of-states description of scanning tunneling spectroscopy 

in supported metal clusters 
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Low-temperature scanning tunneling spectroscopy allows to probe the electronic properties of clus- 
ters at surfaces with unprecedented accuracy. By means of quantum transport theory, using realistic 
tunneling tips, we obtain conductance curves which considerably deviate from the cluster's density 
of states. Our study explains the remarkably small number of peaks in the conductance spectra ob- 
served in recent experiments. We demonstrate that the unambiguous characterization of the states 
on the supported clusters can be achieved with energy-resolved images, obtained from a theoretical 
analysis which mimics the experimental imaging procedure. 

PACS numbers: 73.22.-f, 73.63.-b, 36.40.Cg 



The investigation of the properties of atomic clusters 
is of great fundamental interest since this class of sys- 
tems allows to bridge the gap between single atoms on 
one side and bulk material on the other side [![. This 
applies to both free clusters in the gas phase and clus- 
ters deposited at surfaces. Compared to the gas phase, 
where characteristic effects like magic numbers and size 
dependent bond-character changes are well studied, the 
understanding of supported clusters is much less devel- 
oped due to the the numerous technical difficulties the 
problem poses to both experiment and theory. How- 
ever, there is a growing interest in supported clusters, 
mainly for two reasons: firstly, by means of a scanning 
tunneling microscope (STM) one is now able to perform 
single-cluster studies 2], much more difficult to achieve 
in the gas phase because of the low beam intensity of 
mass-selected clusters. Additionally, scanning tunneling 
spectroscopy (STS) provides valuable information on the 
electronic properties, also allowing spatial resolution, in 
contrast to the non-local character of photoelectron spec- 
troscopy 0. Secondly, supported clusters promise to be 
efficient catalysts as oxidation experiments with small 
gold clusters indicate 0,0. STM/STS techniques repre- 
sent unique tools to analyze supported clusters and thus 
to shed light on the oxidation mechanism. 

STS measurements of metallic clusters (platinum [8|, 
silver 0, and gold dUEl) 011 graphite surfaces were per- 
formed in recent years in order to verify the effect of the 
substrate on the electronic structure of the clusters. The 
STS spectra showed peaks which were unambiguously 
shown to be produced by cluster states. Also, experi- 
ments ruled out that for metal clusters at graphite the ob- 
served peaks are due to Coulomb blockade effects [t| [ll[ . 
Surprisingly, the number of peaks obtained for the con- 
sidered bias- voltage range was considerably smaller than 
the number of electronic states expected for the cluster 
within the corresponding energy interval [1, [t| [13, fill ]. 
These facts clearly indicate that there are states of the 
clusters which do not contribute to the conductance, and 



that the simple interpretation of the STS spectra only in 
terms of the density of states (DOS) [l2j is not sufficient. 
Moreover, no systematic trend of the peak distribution 
as a function of the cluster size could be singled out. Be- 
sides, repeated measurements of the same cluster showed 
variations in the spectra [TlT ] . The blurry picture emerg- 
ing from these observations calls for a theoretical analysis 
in order to evaluate the dependence of the conductance 
spectra on the details of the STS measurement. 

In this Letter we demonstrate that the description of 
STS spectra only in terms of the DOS breaks down in the 
case of metallic clusters on graphite surfaces. Further- 
more, we show that the interpretation and understand- 
ing of STS spectra requires the careful calculation of the 
electron transport through the cluster to the infinite sup- 
porting surface, taking into account a realistic shape for 
the STM tip. Our approach, based on non-equilibrium 
Green's functions (NEGF), treats cluster and surface on 
an equal microscopic footing. It allows us to determine 
the differences between the DOS and the conductance 
for different cluster sizes, tip shapes and surface-tip dis- 
tances. We will show that there is a crucial dependence 
of the STS spectra on the lateral position of the tip which 
should always be kept in mind when discussing experi- 
mental spectra. As we report in our simulations, a reli- 
able characterization of the cluster states should be done 
by means of energy resolve imaging, therefore combining 
spectral information with spatial resolution. 

We performed quantum transport calculations for sil- 
ver clusters having up to 233 atoms on a perfect graphite 
surface. We consider the STM scenario in which the tip 
of the STM is at a distance of a few A over the cluster 
and calculate the tunneling current / between tip elec- 
trode and cluster; all the components are treated at the 
atomistic level. For the cluster we assume the lattice con- 
stant of the Ag fee bulk crystal structure and use a Wulff 
construction [13| to define the facets [14| . By means of 
an effective surface energy for silver and graphite we can 
vary the width-to-height ratio of the clusters to match the 
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experimental value @ . The lateral position of the cluster 
on the surface was found to be of minor importance; the 
cluster-surface distance was set to 2.5 A, consistent with 
ab initio calculations for Ag monomers and dimers on 
graphite [15j. The platinum electrode above the cluster 
is either just a flat ideal surface or a sharp pyramidal Pt 
tip mounted on the electrode surface. The flat electrode 
is used for reference calculations; the sharp tip, having a 
pyramidal structure of three layers with 1-3-6 atoms, can 
be considered as a model of a real STM tip [r| . 

The electronic structure and the transport properties 
of the whole system, i. e. graphite surface, silver cluster 
and platinum electrode/tip, are calculated by means of 
an NEGF approach [53, El. The DOS and the current I 
at a given bias V are calculated from the Green's func- 
tions of the cluster; they contain the interaction with the 
surface and the electrode/tip via tunneling self-energies 
[13]. We employ a self-consistent tight-binding model 
which is parameterized from density-functional calcula- 
tions; atomic charge fluctuations (transfer and polariza- 
tion) are taken into account [l9j ]. The procedure for cal- 
culating equilibrium properties, like the charge transfer 
to or from the cluster, and transport properties, like the 
current through the cluster, is described in Ref. [20j . 

Figure Q] presents a comparison of the DOS with the 
conductance dl/dV for Agn 2 , calculated for both tip 
shapes described above. Two different regions can be 
identified in the DOS of the supported cluster, based on 
the mean level spacing: a rich structure with many peaks 
and level distances of AE w 0.1 eV is visible above the 
Fermi energy (E > 0), while for the occupied states (E < 
0) we have a larger spacing, AE w 0.5 eV. The latter 
value is similar to the energy level spacing measured with 
photoelectron spectroscopy for Ag clusters of comparable 
size, cf. Fig. 2 of Ref. Q. One should mention that, in the 
case of supported silver clusters it is hardly possible to 
assign certain levels to the peaks as it is done for free 
clusters p| or supported fullerenes [2(J. By calculating 
the DOS of the clusters placed at different heights over 
the surface, we observed that, as the clusters approach 
the surface, the levels are not only broadened but are also 
shifted to higher or lower energies in a complicated way. 

Next, we compare the DOS to the conductance pat- 
terns for the flat and the sharp tips, taken as represen- 
tatives of two limiting cases. Panel (b) of Fig.Q] shows 
dl/dV for the flat electrode. The relative peak height 
changes compared to the DOS; some peaks are reduced, 
for example the peak A at Ea— — 1.18 eV, while oth- 
ers are enhanced, for example at E-q= — 0.54 eV. Despite 
these variations, the overall picture is largely preserved 
and we conclude that the conductance measured with the 
artificially flat electrode would give a reasonable picture 
of the DOS of the cluster. The situation changes dras- 
tically for the sharp tip, as shown in panel (c) of Fig.Q] 
Due to the reduced number of atoms of the probe partic- 
ipating to the tunneling process, the conductance drops 
by more than an order of magnitude (in fact we observed 
that for the flat electrode the average conductance scales 



directly with the nuber of atoms in the top facet of the 
cluster) . More remarkably, a few peaks completely disap- 
pear, for example the one labeled A, while many others 
are strongly suppressed. The discrepancy between the 
spectra obtained with the two kinds of electrodes can be 
understood as follows. While the flat electrode signal re- 
ceives contribution from all atoms on the top facet, repre- 
senting the average DOS, the sharp tip probes the local 
DOS in a small region, and the resulting conductance 
will be affected by the spatial dependence of the specific 
state (see below). Independently of the detailed struc- 
ture of the probe in the experiments, we expect that the 
measured conductances should be related to our sharp- 
tip results. Any atom sticking out of the tip will carry 
almost all current in the tunneling regime. 

We have extended the study to other cluster sizes and 
present six of them in Fig. [5] In each of the graphs we 
compare the sharp-tip conductance with the DOS of the 
supported cluster. As in Fig.Q]we placed the tip centrally 
over the cluster. Additionally, the conductance curves 
were convoluted with a Gaussian function of 100 meV full 
width at half maximum, in order to reproduce the accu- 
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FIG. 1: (Color online) Density of states (panel a) of an Agn2 
cluster at a graphite surface compared to the conductance 
spectra obtained with flat (b) and sharp (c) platinum elec- 
trodes as a function of energy, being E — the Fermi energy. 
Go = 2e 2 /h is the quantum unit of conductance. 
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FIG. 2: (Color online) Conductance (thick blue lines, left axis) through silver clusters of various sizes on a graphite surface; 
the sharp Pt conductance curves are convoluted with a Gaussian having 100 meV full width at half maximum. For comparison 
the density of states (thin lines, right axis) is shown as well. Go = 2e 2 /h is the quantum unit of conductance. 



racy of experimental low-temperature measurements [ll| . 
Even if the detailed shape of the curves differs from one 
cluster size to the next one, the general message is the 
same: the conductance is much less structured than the 
corresponding DOS. Thus we conclude that the experi- 
mentally found larger level spacing in supported cluster 
[1, 0, [l(| [HI, if compared to the gas phase [H, is due 
to the local and selective character of the transport mea- 
surement. However, in accordance with the experimental 
observation , positions and widths of the conductance 
peaks in Fig. [2] do not show any trend with the cluster 
size. 

So far we have described situations in which the sharp 
tip was placed over the cluster center. Moving the tip lat- 
erally from the central position may change the observed 
conductance spectra substantially. Energy resolved maps 
are obtained by projecting the local conductance on a 
constant current isosurface at particular energies Ex- 
The constraint of a constant current isosurface is nec- 
essary to compare conductance values measured at dif- 
ferent positions over the sample. This point represents 
the main difficulty in the calculations, since tip trajectory 
invariably describes a complicated surface. 

In Fig. [3]we report such maps for the states labeled A- 
D in Fig.[TJ Figure [3] shows for each of the chosen energies 
the constant-density isosurface (upper row) and the con- 



ductance map (lower row). The former has been obtained 
by integrating the lesser Green's function over a small en- 
ergy interval of 20 meV centered around the energies Ea, 
Eb, Eq, and Ed in order to get the charge density in 
terms of the local atomic orbitals. The constant-density 
isosurfaces show that all four states are delocalized over 
the whole clusters. However, their shapes differ, in par- 
ticular at the upper facet. 

A much clearer picture of these different structures 
emerges from the conductance maps, shown in the lower 
row of Fig. [3] These maps display as a contour plot 
the conductance dI/dV(x, y, z, Ex)\i=i on a constant- 
current isosurface z = z(x,y) determined by the implicit 

equation I(z; x, y) = Iq. In the calculations, the tip was 
moved on a three-dimensional finite-element mesh with 
a grid spacing of 0.5 A. At each grid point the current 
and the conductance were computed, and the constant- 
current isosurface was obtained by interpolation. Having 
this three-dimensional topographic image of the cluster 
we render conductance maps by showing the calculated 
dl/dV values with a colored contour plot on it. The 
obtained patterns are dominated by the three-fold sym- 
metry of the clusters, which is, however, disturbed by the 
coupling to the supporting surface. The central minima 
(blue spots in online version, black spots in print ver- 
sion) of the patterns A and C clearly reveal the reason 
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FIG. 3: (Color online) Constant-density isosurfaces (upper row) and conductance maps (lower row) for four selected states of 
Agn2. The four peaks are marked by red lines in Fig.[T] their energies are Ea=— 1.18 eV, Eb=— 0.54eV, -Ec=+0.42 eV and 
Er>—+l-20eV, respectively. The conductance maps are presented as contour plots on constant-current isosurfaces (see text); 
red (light grey) = high conductance, blue (dark grey) = low conductance. The side length of the squares is 14.5 A. 



for the disappearance of the corresponding peaks in the 
measurement with the tip over the centre, cf. Fig.[TJ Fur- 
thermore, it becomes clear that placing the tip a few A 
off the center leads to considerable changes in the con- 
ductance of these peaks. This may explain the observed 
discrepancies in the spectra obtained from measurements 
performed on the same cluster but at different times pd| . 

Summarizing, we have shown that recent data from 
STS measurement of metallic clusters at surfaces can be 
understood only in terms of non-equilibrium transport 
calculations. The calculated current and conductance, 
respectively, depend crucially on the tip shape. A sharp 



tip is very selective, which explains the low number of 
states seen in the experiment compared to the richly 
structured DOS of the cluster. This selectivity depends 
strongly on the tip position, which leads to the require- 
ment of a more complete experimental characterization 
based on energy-resolved imaging of the electronic states 
on the cluster [2l| . 

We acknowledge financial support by the Deutsche 
Forschungsgemeinschaft through the priority program 
SPP 1153 "Clusters in contact with surfaces: Electronic 
structure and magnetism" . 
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